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Persistent scaling behavior of magnetization in layered high T c superconductors with short-range 
columnar defects is explained within the Ginzburg Landau theory. In the weak field region, the scal- 
ing function differs from that of a clean sample and both the critical and crossing temperatures are 
renormalized due to defects. In the strong field region, defects are effectively suppressed and scaling 
function, as well as critical and crossing temperatures are the same as in a clean superconductor. 
This picture is consistent with recent experimental results 



Layered high-temperature superconducting (HTSC) 
materials, such as Bi2Sr2CaCu20s+<5 and Bi2Sr2Ca2Cu3- 
Oio, are known to exhibit 2D scaling magnetic proper- 
ties around the mean field transition line H C 2(T). 
It is manifested by inspecting the magnetization Mq as 
a function of temperature T (measured in energy units) 
and the (external) magnetic field H B : 



^L=M Q (T,H) = -(V~* 



(1) 



where s is an effective interlayer spacing, $o is the flux 
quantum, x = AH' c2 [T - T c2 (H)]/VTH is the scaling 
variable, H' c2 = —dH C 2(T) / cIT\t=t c0 , and T c q is the zero 
field critical temperature. For a superconductor with 
Ginzburg-Landau (GL) parameter k and Abrikosov ge- 
ometric factor (3a the constant A = y^s^o/p, where 
p = 16ttk 2 /3a- The form of the scaling function Q implies 
the existence of a crossing point: at some temperature 
Tq = T c0 (l + \/{2A 2 H' c2 ))- 1 , the sample magnetization 
is independent on H, Af * = Mq(Tq, H) = -T */(s$ ). 

Recently, the influence of linear defects (columnar de- 
fects, artificial holes etc.) on the magnetic properties of 
superconductors has been studied experimentally J^-j?]] 
and theoretically |8| [L2| . In particular, experiments by 
van der Beek et. al. J7[ showed that in HTSC with colum- 
nar defects the reversible magnetization of the sample 
is drastically affected, and that there are now two scal- 
ing regimes pertaining to relatively weak H < H<& and 
strong H > iJ$ magnetic fields (here the matching field 
H<j> = rid&o is proportional to the 2D density of defects 
rid). These two scaling regimes correspond to two dif- 
ferent critical temperatures (used in Ref. (7) as fitting 
parameters) and crossing points. 

In this Letter we propose an explanation of these re- 
sults. Let us commence by presenting some intuitive ar- 
guments. Consider the quantity c = H$/H which, in a 
macroscopic sample, is the number of defects divided by 
the number of vortices. The magnetic field then serves 
as a control parameter for tuning the effective concen- 
tration c of defects. In the weak field region, c is large, 
each vortex is affected by a force emanating from many 



defects, and the fluctuations of this force play the main 
role. Short-range defects could be taken into account 
perturbatively. In first order they retain the same form 
of scaling function as that of a clean sample but renor- 
malize the critical temperature T c . Second order correc- 
tions indeed destroy the scaling behavior but in the vicin- 
ity of the crossing temperature scaling is approximately 
maintained. In the strong field region, c is small, and 
the standard concentration expansion jl3| can be used. 
Here, even the first order correction (with respect to 
small concentration) destroys the scaling behavior. How- 
ever, a strong field effectively suppresses the defects, thus 
restoring the scaling behavior of a clean superconductor 
with the initial critical temperature T c q. Identifying the 
two fitting temperatures of Ref. Q with the renormal- 
ized critical temperature T c and the initial one T c q re- 
spectively, one finds for the dimensionless defect strength 
6i = 0.49, well inside its allowed range < B\ < 1. This 
indicates a full consistence between the description con- 
structed below and the experimental results of Ref. . 

Our quantitative discussion employs an approach pro- 
posed and successfully used for arbitrary fields in clean 
superconductors |^|,|],[l4j and for very low fields in dis- 
ordered superconductors JhJ. Here we use it for dis- 
ordered superconductors in much higher fields. Con- 
sider an irradiated thin superconducting film (or one 
layer in a layered superconductor) with area S subject 
to perpendicular magnetic field (thus parallel to the de- 
fects). The effective interlayer separation s is assumed 
to be much larger than the effective superconducting 
coherence length £,(H, T) in the magnetic field direc- 
tion but much smaller than the magnetic penetration 
depth. Then the problem becomes effectively two di- 
mensional juj. Columnar defects can be described as 
a local reduction of the critical temperature 5T c (r) = 
Tea tj ex P(~ ( r — r j) 2 /2£ 2 )- Here r is a two dimen- 
sional vector in the film plane, L is the defect radius, 
and the positions Tj of defects are uniformly and inde- 
pendently distributed over the film plane with density rid- 
The value of rid is assumed to be moderate so that for the 
pertinent region of temperature the matching field H§, 
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is always much smaller than H C 2(T). The dimensionless 
amplitudes of defects tj < 1 are also independent ran- 
dom quantities distributed with some probability density. 
The thermodynamic properties of a type-II superconduc- 
tor with k»1 containing N v vortices arc described by 
its partition function 



Zcx J V{9}exp(-N v g[9]), 



(2) 



where \& is the corresponding order parameter. The di- 
mensionless GL free energy g[^f] of an irradiated super- 
conductor is given by an expression 



g = x\^ + (4p A )-^ + T\^ 



(3) 



where bar denotes averaging over the sample area. The 
scaling variable x and the local temperature r(r) are de- 
fined below for each region of the magnetic field. 

Following ITS], we replace in Eq.(||) \^>(r)\ by 

Pa ^^(r)! 2 ^ * where Pa ~ 1.16 is the Abrikosov fac- 
tor for a triangular lattice. This replacement is based on 
the assumption that the distribution of vortices is almost 
uniform in both regions of the magnetic field considered 
here. It is supported by noticing a remarkable difference 
between the number of vortices and the number of defects 
in both regions of fields |15|| . This substitution, together 
with the simplest version of the Hubbard-Stratonovich 
transformation (introduction of an additional integration 
over some auxiliary field 7) turns the problem to be an 
exactly solvable one ]l0| . Then project the order param- 
eter on the lowest Landau level (LLL) subspace, 



*(r) 



£ 

m=0 



(4) 



where L m (r) are normalized LLL cigenfunctions with or- 
bital momentum m. As was recently demonstrated [ OJ , 
the LLL approximation works quite well even down to 
H > H C 2(T)/13. After integration over the expansion 
coefficients C m the partition function (0) reads, 



Z oc J exp{—N v L(j,x)}dj, 

where 

£(7,x) = -7 2 + AT-Hrln (x + 7)I + f 
and f is a random matrix with elements: 



/ L* m {v) T {v)L n {v)d 2 v. 
Js 



(5) 



(6) 



(7) 



The contour T in Eq.(||) is parallel to the imaginary 
axis and stretches from 7* — ioo to 7* + too. To as- 
sure convergence of the integrals over the coefficients 



{C m } the real constant 7* should satisfy the inequality 
7* + x + minT„ > 0, where r« is the n-th eigenvalue of 
the matrix r m „. 

In the thermodynamic limit S — » 00 with = N v /S 
fixed, the partition function (|J) could be calculated in a 
saddle point approximation. This results in the following 
form for the magnetization 



-.M{T, H) = {N V Z) 1 dZ/dx = -2rf(x), (8) 



where j(x) is the solution of the saddle point equation 
dC( r ),x)/d^/ = 0. For a clean superconductor (f = 0) 
one gets two possible saddle points but only one of them 



7o (x) = -(Vx 2 + 2-x) 



(9) 



can be reached by an allowed deformation of the contour 
T. Substitution of Eq.(g) into (||) yields the magneti- 
zation Mq(T,H) of a clean sample (|l|) obtained in Ref. 
0. Note that — 270(2;) serves as the appropriate scaling 
function. To study the disordered case, we consider sep- 
arately two regions of the magnetic field. 

In the weak field region H < H$ we, from the on- 
set, take into account the renormalization of the critical 
temperature caused by defects. As a result, the scaling 
variable x is defined in the same way as for a clean super- 
conductor albeit with renormalized critical temperature 
T c = T c0 - ST Cl where 



ST C =< ST c {r) >= 2TT0 1 n d L 2 T l 



cO 



(10) 



and 6 n =< t n > where here and below < .. > implies 
ensemble average. The function r(r) in this field region, 
defined as r(r) = (5T c (r) - ST c )AH' c2 /VTH, represents 
temperature fluctuations caused by short-range defects. 
They are small and can be accounted for perturbatively. 
Then, in the thermodynamic limit, the last term on the 
r.h.s. of Eq. (|6j) has an explicit self-averaged structure 
A r I 7 1 tr(...) and can be replaced by its average. This pro- 
cedure modifies the saddle point equation and therefore 
results in a modified magnetization 

M(T,H) = M (T,H) fl + £ (T)4==,") , (H) 



V?T2' 



where 



clT) =(^y ^*w.L* (12) 

Note that the parameter e(T) is proportional to the 
fourth power of the defect radius L thus justifying the 
perturbation approach for short-range defects. 

In the zcroth approximation with respect to e(T) the 
magnetization ( |Tl| ) has exactly the same form, as for a 
clean sample (Q) thus retaining both the scaling property 
and the existence of a crossing point. However, due to 
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renormalization of the critical temperature, the crossing 
temperature T* = T * - 5T* differs from its value T * in 
a clean sample: ST* = ST C (1 + {2A 2 H' c2 )- 1 )~ 1 . In the 
next order, scaling is virtually destroyed, since the cor- 
rection term (within the parenthesis in Eq.(|ll])) depends 
not only on the scaling variable x but also on temper- 
ature. Yet, in a sufficiently narrow region around some 
temperature T, the deviation from scaling is negligibly 
small, but the scaling function itself is modified to be 
—2 r y(x,T). At temperature T* the magnetization reads 



M{T*,H) = M (T*) 1 + e(T*) 



2H* 



H + H* 



(13) 



where if* = H c2 (T*) = T*/(2A 2 ). Therefore if the field 
is weak enough, if <C if* , then the crossing point is re- 
stored, T* serves as a true crossing temperature and the 
magnetization at the crossing temperature differs from 
its unperturbed form —270(2;) merely by a multiplicative 
constant 1 + 2e(T*). 

When the magnetic field increases, the approach used 
above becomes inapplicable. Firstly, it fails in the vicin- 
ity of the matching field where the Abrikosov factor be- 
comes very sensitive to the details of defect configura- 
tion. Secondly, higher order terms in the perturbation 
expansion for the saddle point equation (which are omit- 
ted), grow with magnetic field. Fortunately, we have 
here a new small parameter, that is, the dimensionless 
concentration c of defects. It is then natural to use 
the concentration expansion. In such a case there is 
no sense in renormalizing the critical temperature, and 
the dimensionless temperature r(r) is now defined as 
r(r) = 6T c (r)AH' c2 /VTH. 

As mentioned above, the second term in the r.h.s. of 
Eq.(^j) is self-averaging and can be calculated using the 
limiting form of the density of states p(r) of the matrix 
(Q), which, for short-range defects in linear approxima- 
tion with respect to c, reads 



P (t) = (1-c)*(t) + 



(14) 



where A = 2nL 2 T rM AH' c2 v r H(<P VT)- 1 and p(t) is prob- 
ability distribution of the dimensionless temperature tj . 
Indeed, the matrix r mn is nothing but the Hamiltonian 
of a particle with charge 2e in a 2D system subject to a 
perpendicular magnetic field and containing short-range 
defects (projected on the LLL). The first and second 
terms in Eq. ( |l4| ) correspond, respectively, to those states 
whose energy is stuck to the LLL (despite the presence of 
zero-range defects (see e.g. and those states whose 

energies are lifted from the LLL by these defects. For 
sufficiently narrow distribution p(t), the corresponding 
saddle-point equation leads to the magnetization 



M = M 1 



c\9 1 



(1 + 2Xe ll0 (x))Vx 1 



(15) 



were Mq(T, if) is given by Eq.(Q) with an initial critical 
temperature T c q. 

Rigorously speaking, scaling is destroyed since both 
the concentration c and the shifted eigenvalue 9\X de- 
pend explicitly on if and T. However, at strong field 
the correction term in Eq. ( [l5| ) becomes negligibly small. 
This implies a restoration of the crossing point. Indeed, 
at temperature T * the magnetization M* = M(T*,H) 
assumes the form 



M* = MX 1 - 



1 



ifj> 



1 + 77 if + if * 



(16) 



with rf l = 2TrL 2 H' c2 T cQ 9i/§ . Therefore in the entire 
strong field region ff$ <C if <C if* the crossing tempera- 
ture coincides with its initial value T * and the magneti- 
zation in the crossing point practically coincides with its 
value Mq in a clean superconductor. 

Let us now discuss the limits of applicability of our 
results and their relation to the experiment of Ref. 
0. Note that the first two moments 0\ t2 of the ran- 
dom dimensionless temperature t satisfy the inequality 
< Of < 62 < 1. In the pertinent region of fields 

0.2 5T, a typical defect radius L ~ 3.5nm is at 
least one order of magnitude smaller than the magnetic 
length, hence the defects can definitely be taken as short- 
range ones. In the weak field region, the important small 
parameters are then e(T) (which enters the magnetiza- 
tion (p"3|)) and e(T)/(x + jq(x)) 2 (which enters the sad- 
dle point equation). Using parameters from the exper- 
imental setup s = 1.5nm, ks[J<oH' c2 = 1.15TK -1 , 
k = 100, n d = 5 x 10 10 cm- 2 , /j ff$ = IT, T = 75H-85K, 
T* = 78. 9K, wc find from Eq.@ e(T*) = 0.56> 2 and 
e{T*)/(x* +~/ {x*)) 2 w 0.25 (the latter figure is obtained 
for fi Q H — 0.2T). For quite plausible value 9 2 = 0.5 one 
then finds s(T*) = 0.56> 2 = 0.25. The condition of con- 
vergence of the integral over the expansion coefficients 
{Cm] can be written as if > Q.2bH^,9 2 /0 2 and even in 
the worst case Q\ = 82 it reads /Joff « 0.25T. Finally, one 
has /Joff * ~ 6.4T and applicability of the LLL projection 
requires [IqH > 0.5T. The weak field region of Ref. 
corresponds to /j,qH — 0.2 -j- 0.02T. Thus, in the weak 
field region, the condition for applicability of the LLL 
projection is slightly violated, but the deviation is not 
dramatic. In the strong field region we find 77 ps 2.9 and 
therefore the correction term in parenthesis of equation 
( |l6| ) is less than three percents. Hence, in this region our 
assumptions are fully satisfied. 

Using the same set of parameters we display in figure 

1 the quantity Mj y/TH as a function of the scaling vari- 
able for weak field (inset) and strong field (main part). 
We used here the maximal value d\ = 1. In the strong 
field region, the deviation form clean sample scaling be- 
havior is negligibly small for all three values of strong 
magnetic field. 
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FIG. 1. The quantity M/V HT, as a function of x. Dashed, 
dotted and dot-dashed lines correspond to fioH = 3T, 4T 
and 5T (main figure) and correspond to fioH — 0.02T, 0.1 
and 0.2T (insertion). The solid line corresponds to the 
clean-sample scaling function (strong field region only). 

In the weak field region, the scaling functions for three 
different fields can hardly be distinguished. This means 
that scaling is undoubtedly valid in a vicinity of the cross- 
ing temperature. At the same time the scaling function 
differs from its form in a clean sample (|l|) by a multi- 
plicative constant (see the parenthesis in Eq.(JTl|)). Note 
that scaling in the weak field region (which was exper- 
imentally established) is less pronounced than that in 
the strong field region. Apparently, the reason is that 
the experimental data are fitted to account for the clean 
sample scaling function. Nevertheless if we identify the 
fitted temperature 82. 6K (found in Ref. pj] in the weak 
field region) with the renormalized critical temperature 
T c = T c q — 5T C , and the fitted critical temperature 84. 2K 
in the strong field region with T c o, then, even within such 
a rough approximation, we obtain 6\ « 0.5. Recalling 
that 0\ should be positive and less than unity, the above 
result strongly supports the applicability of our theory 
to the pertinent experiment Q . 

In summary, we calculated the magnetization of an 
irradiated superconductor below the mean-field transi- 
tion line iJ c2 (T), using the approach developed in Refs. 
P,^| JlO|Jl4]| . It was shown that, from a rigorous point of 
view, disordered short-range defects are expected to de- 
stroy the scaling behavior and prevent the existence of 
crossing point in both regions of weak and strong mag- 
netic fields (with respect to matching field H<&). And 
yet, in the framework of the experimental setup [M the 
deviation from scaling behavior appears to be negligibly 
small and crossing points exist in both field regions, in 
complete agreement with the experimental findings. The 
two fitting critical temperatures introduced in Ref. |7]] 
for the strong and weak field regions correspond, in our 
formalism, to the initial and renormalized critical tem- 
peratures. 

The authors would like to thank Z. Tesanovic for help- 
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